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Pli' abstract 

We have used the 76-m Lovell, 94-ni equivalent WSRT and 100-m Effelsberg radio 
5^ . telescopes to investigate the simultaneous single-pulse properties of the radio emitting 

^ ■ magnetar AXP XTE J1810-197 at frequencies of 1.4, 4.8 and 8.35 GHz during May 

I and July 2006. We study the magnetar's pulse-energy distributions which are found 

to be very peculiar as they are changing on time-scales of days and cannot be fit by 
a single statistical model. The magnetar exhibits strong spiky single giant-pulse-like 
subpulses, but they do not fit the definition of the giant pulse or giant micropulse 
phenomena. Measurements of the longitude-resolved modulation index reveal a high 
^vjj I degree of intensity fluctuations on day-to-day time-scales and dramatic changes across 

QQ , pulse phase. We find the frequency evolution of the modulation index values differs 

fT^ ' significantly from what is observed in normal radio pulsars. We find that no regular 

drifting subpulse phenomenon is present at any of the observed frequencies at any 
observing epoch. However, we find a quasi-periodicity of the subpulses present in 
' the majority of the observing sessions. A correlation analysis indicates a relationship 

, between components from different frequencies. We discuss the results of our analysis 

in light of the emission properties of normal radio pulsars and a recently proposed 
model which takes radio emission from magnetars into consideration. 

I Key words: magnetars: general - stars: individual: AXP J1810-197 ~ stars: neutron 

■ - pulsars 



1 INTRODUCTION 

Anomalous X-ray Pulsars (AXPs) and Soft Gamma-ray Re- 
peaters (SGRs) form a class of slowly rotating neutron stars 
with rotational periods ranging from 2 to 12 seconds, called 
magnetars. Magnetars are characterised by properties such 
as: emission over a wide spectrum of wavelengths (7-ray, X- 
ray, optical, near infrared and radio), rapid spin-down and 
a very high magnetic field (higher than 10^* G, which is two 
orders of magnitude higher than for a radio pulsar of the 
same age and three orders of magnitude greater than for 
an average radio pulsar) . These properties are very peculiar 
compared to the normal rotating neutron star behaviour and 
are best underst ood in the context of the ma gnetar model 
first discussed bv lDuncan fc ThompsonI (|l992l ). 
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The AXP XTE J1810-197 wa s serendipitous l y disc ov- 
ered in 2003 as an X-ray pulsar bv llbrahim et all l|2004l ) in 
data taken with the Rossi X-Ray Timing Explorer (RXTE) 
while observing the known Soft Gamma-ray Repeater SGR 
1806-20. A search in the archival XTE data showed that it 
produced an outburst around November 2002 with a per- 
sistent decline of its X-ray fiux since then. The reported 
X-ray pulsar spin period was 5.54 s with a spin-down rate 
~ 10"" ss"^ which, using standard methods, implied a 
magnetic field strength of ~ 3x lO" G ('Ibr ahim et"aLll2004l '). 
Due to this unusual long-term fiux variability (outburst and 
exponential decline) it was classified as the first transient 
magnetar. 

In 2004 a radio source with a strong fiux density of 
4.5±0. 5 mjy was found at t he exact position of XTE J1810- 
197 bv lHalpern eTall (|2005h in data from the VLA MAGPIS 
survey at 1.4 GHz. Archival data were searched for any ear- 
lier detections but only upper limits could be estimated. Ob- 
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servations made bv lCamilo e\~el\ (|2006l ) at the Parkes radio- 
telescope and the Green Bank Telescope discovered the pres- 
ence of the first pulsed radio emission from a magnetar. In a 
following series of papers, Camilo and collaborators investi- 
gate the magnetar's properties including: polarisation, flux 
density and pulse mo rp hology variations in the radio regime 
l|Camilo et al.ll2007bl & ICamilo et allbOOTal) and broadban d 
observations (radio, infrared to X-ray; Camilo et al.ll2007j ). 
The results show large variations in pulse morphology and 
flux with an overall tendency for the flux to decrease during 
mid-2006 observations. In the latter paper the authors re- 
port the detection of single pulses at 88 GHz with the IRAM 
30 m telescope and an infrared counterpart with variable flux 
and is uncorrelated with either X-ray or radio fl uxes. 

In their r ecent papers, iKramer et al.l l|2007h and 
iLazaridis et al.l l|2008l) presented results from the first, full 
polarisation, simultaneous multi-frequency observations of 
XTE J181Q-197. Although the magnetar shares some prop- 
erties observed in ordinary radio pulsars, some distinct dif- 
ferences were also found. The most peculiar is the very 
flat and variable spectrum with an average spectral index 
a — 0.00 ±—0.08, which makes it the brightest neutron star 
emitting at frequencies above 20 GHz. Another difference 
is the variability of the pulse profile. Its extreme changes 
of shape on time-scales of less than a day, makes it quite 
unusual when compared to normal radio pulsars which typ- 
ically need only several hundred individual pulses to obtain 
a stable pulse profile. Moreover the high degree of linear po- 
larisation in the main pulse (~ 90%) and the position angle 
evolution with time provides evidence of a clear difference 
in the magnetar emission mechanism compared to normal 
radio pulsars. 

Thanks to the initially extremely bright nature of XTE 
J1810-197 we are able to detect its single pulses, by analysis 
of which, we will try to understand the long-term pulse mor- 
phology and if it shows any correspondence with the single- 
pulse properties of normal radio pulsars. We report simul- 
taneous multi-frequency single-pulse observations conducted 
at frequencies of 1.4, 4.9 and 8.35 GHz during observing ses- 
sions in May and July 2006. We present the results of our 
search for subpulse modulation of XTE J1810-197. We also 
investigate the pulse-energy distributions at all frequencies 
and epochs and the pulse-to-pulse correlation properties be- 
tween the simultaneous single pulse observations. We then 
consider these properties in light of the emission character- 
istics of normal and giant pulses from ordinary pulsars. 



2 OBSERVATIONS 

The observations were made using the 94-m equivalent 
Westerbork Synthesis Radio Telescope (WSRT) in the 
Netherlands, the 76-m Lovell radio telescope at Jodrell Bank 
Observatory of the University of Manchester, UK and the 
100-m radio telescope of the Max-Planck Institute for Ra- 
dioastronomy (MPIfR) at Effelsberg, Germany, which are 
working in the European Pulsar Network (EPN) collabora- 
tion. The detailed discussion of the o bserving systems an d 
calibration procedures is described bv lKramer et al.l (|2007h . 
while Table [T] summarises the details of the observing ses- 
sions used in this paper. 
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Figure 1. An example of an average pulse profile of XTE J1810— 
197 from the observation made with the Effelsberg radio telescope 
during session 3. This plot shows the component naming conven- 
tion used in this paper. 



3 ANALYSIS AND RESULTS 

Observations of individual pulses from XTE J1810-197 al- 
lowed us to perform a variety of analysis techniques which 
we describe below. Before performing the analysis, all data 
sets have been corrected for any RFI present. In the case of 
the Lovell and WSRT observations re-binning was applied 
to increase the signal-to-noise ratio and to match the 5.4 
ms time resolution of the Effelsberg data set. If more than 
one dataset was present per telescope per session, they were 
aligned in longitude and appended together to improve the 
statistics. This was done only for pulse-energy distributions, 
because any interruption (i. e. missing pulses) between con- 
secutive data sets affects fiuctuation and correlation analysis 
results. For the correlation analysis we aligned the data in 
phase and pulse number. Starting from session 2 we iden- 
tify two emission regions in the average pulse profile, la- 
belled as main pulse (MP), and in t erpuls e (IP) according 
to the classification of lKramer et al.l l|2007l ). The MP region 
can be divided into at least three separate longitude regions 
containing separate components as marked in the panels in 
Fig. [1] The average profile from the 1.4 GHz Lovell data 
in session 1 and 2 lacks the third component visible in the 
4.9 GHz WSRT and 8.35 GHz Effelsberg data sets. The IP 
is not present in the first session and is only visible in the 
Effelsberg data in the second session. It is however visible in 
all three data sets in the third session. For the full list when 
the components are present see Table [1] 



3.1 The shape and stability of the pulse profile 

Fig.[2]presents an example of a sequence of pulses from XTE 
J1810-197 observed at a frequency of 8.35 GHz during ses- 
sion 1. The plotted pulse- longitude range corresponds to the 
whole range of MP as can be seen in Fig. [T] One can easily 
see that the subpulses are much narrower than the width of 
the average profile and appear at the longitude ranges corre- 
sponding to different components (see for instance the sub- 
pulses around pulse number 50). Also, it is worth noticing. 
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Table 1. Summary of observing sessions. 



Session Date Telescope Frequency BW Total number Component 





MJD/dd.mm.yy 




[GHz] 


[MHz] 


of pulses 




presence 




1 


53886/31.05.06 


Lovell 


1.418 


32 


2101 


Ml 


M2 - 








WSRT 


4.901 


80 


2151 


Ml 


M2 M3 








Effelsberg 


8.350 


1000 


972 


Ml 


M2 M3 




2 


53926/10.07.06 


Lovell 


1.418 


32 


1947 


Ml 


M2 - 








Effelsberg 


8.350 


1000 


2583 


Ml 


M2 M3 


IP 
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IP 
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Figure 2. An example of sequence of successive single pulses from 
the MP phase of XTE J1810-197 from the observation made at 
8.35 GHz with the Effelsberg radio telescope during session 1. 

that the subpulses associated with Ml are stronger than the 
subpulses from the remaining MP components, at any ob- 
served frequency and epochs, except the 8.35 GHz Effelsberg 
data set in session 3 where M2 becomes the strongest com- 
ponent in the MP longitude range (see Fig. The strong 
and spiky subpulses tend to be separated by ~ 4 deg (~ 61 
ms) throughout the data. 

3.2 Fluctuation analysis 

It was first shown by brake fc Craftl { 19681 ) that for some 
pulsars, subpulses exhibit 'drifting' across the pulse win- 
dow in an orderly fashion. This phenomenon forms 'drift 



bands' as can be seen in the 'pulse-stack' formed when one 
takes successive pulses and plots them on top of one other. 
The techniques used in this paper to investigate this phe- 
nomenon involve the Fourier transform and are kn own as 
the lo ngitude-resolved fluctuation spectrum (LRFS; iBackej 
197Cf) and two-dimensiona l fluctuation spectrum (2DFS; 
Edwards fc Stapperd |2002| ) . The procedure for the fluctu- 
ation analysis of our data sets is identical to those presented 
bylWeltevredc et al. (2007), so we will present only a basic 
summary here. 

The first step in the fluctuation analysis is to form an 
average profile from a pulse-stack by vertically integrating 
each phase bin within the same longitude along consecutive 
pulses. Fig. [3] shows the results from analysis at frequencies 
of 1.4, 4.9 and 8.35 GHz from session 1, where the aver- 
age profile is drawn with the solid line in the top panel. 
The abscissa denotes the pulse phase in degrees. The top 
panel also presents the longitude-resolved modulation index 
(LRMI; solid line with error bars). The longitude-resolved 
modulation index is a basic method to estimate the presence 
of subpulse modulation. We only show LRMI values when 
the error in the LRMI is less than 0.5. 

To investigate whether this modulation is systematic or 
lacks organisation, the 2DFS has to be calculated. This is 
done by dividing the pulse-stack into blocks of 512 pulse^ 
and applying the Discrete Fourier Transform (DFT) along 
lines with different slopes in the pulse-stack. After averaging 
the spectra from each transformed block the final spectra 
is produced (Fig. |31 panels below average profiles). If the 
pulsar exhibits subpulse drifting this will be visible in its 
spectra as a region, so-called feature, of enhanced power 
in the greyscale. The vertical separation of possible drift 
bands is denoted by P3 which is expressed in pulsar periods, 
Po- The information about the horizontal separation of the 
possible drift bands, and thus whether the subpulses are 
drifting over a certain longitude range, is denoted by P2 and 
is expressed in longitude units. A positive or negative P2 
value of a feature in the 2DFS denotes the separation of 
the subpulses, indicates that the subpulses appear later or 
earlier in successive pulses respectively. 

We present the average profiles, LRMI and 2DFS from 
all telescopes from session 1 in Fig. [3] Results from our fiuc- 

^ Wherever the number of pulses was sufficient, otherwise shorter 
transforms was used. 
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Figure 3. Fluctuation analysis results shown for frequencies of 1.4 GHz {left), 4.9 GHz (middle) and 8.35 GHz {right) from session 1 
(MJD 53886). The upper panels show the integrated pulse profile (solid line) with the peak amplitude normalised to 1, and the longitude- 
resolved modulation index, LRMI (solid line with error bars). The lower panels show the 2DFS where the ordinate of the resulting spectra 
are given in cycles per period (cpp) which corresponds to Pq /P3 (Pq is the pulsar period and P3 denotes the vertical separation of possible 
drift bands) and the abscissa is also in cpp but now it corresponds to Pq /P2 (P2 denotes the horizontal separation of possible drift bands 
in longitude units). The greyscale intensity of 2DFS corresponds to the spectral power. The presence of a significant spectral feature 
with a value of P3 and a positive or negative value of P2, indicates that the subpulses appear with a preferred periodicity of P3 pulsar 
periods and horizontal separation given in longitude units. The side panels correspond to horizontally (left panel) and vertically (bottom 
panel) integrated spectrum. 



tuation analysis show no visible spectral features indicating 
that there is no regular drift patterns in any of our observa- 
tions. However, one can see that there are spectral features 
in the 2DFS plots indicating subpulse modulation not asso- 
ciated with drifting subpulses, which will be discussed below. 
The LRMI in the 8.35 GHz Effelsberg data set shows a sig- 
nificant drop in the middle of Ml, while towards its trailing 
edge the LRMI values increase. Inspecting the pulse stack 
indicates that this drop is caused by the frequent occurrence 
of many strong and narrow subpulses with similar intensi- 
ties at that pulse phase. Similar behaviour appears in all the 
data sets in this session, but it is less visible in the 1.4 GHz 
Lovell and the 4.9 GHz WSRT data sets. In both the WSRT 
and the Effelsberg data sets the LRMI values are larger than 
the Lovell data set by a factor of two. The vertically inte- 
grated 2DFS (lower panels in Fig.|3)) show similarity in shape 
across all the frequencies but change in intensity. There are 
two bumps and a strong peak in the middle of the vertically 
integrated 2DFS. We interpret the bumps as apparent sub- 
pulse separation, i.e. subpulses tend to be equally separated 
in the pulse profile throughout the observation. The peak 
is not related to the magnetar behaviour but indicates a 
baseline variations due to effects extrinsic to the source and 
instrumental effects due to the weather. This peak is also vis- 
ible in the vertically integrated 2DFS from the WSRT data 
set, but it is weaker. In the Effelsberg data set the peak is 
no longer visible, while the bumps increase in intensity. 

In session 2, the 2DFS plots show no visible spectral 
features indicating there is no preferred periodicity. The av- 
erage profiles indicate small changes from session 1 to session 
2 indicating that the subpulse separation is now less regu- 
lar. The 2DFS plots made from the 8.35 GHz Effelsberg 
data sets for both the MP and IP show a strange wave- 
like pattern (Fig. |4l middle panel) . Further analysis showed 



that this rapid pattern is due to interference and may be 
caused by the convolution of the interfering signal from the 
alternating stray magnetic fields of the fans used to cool the 
telescope receivers with the actual on-pulse magnetar sig- 
nal. The pattern is stronger in the IP 2DFS plots due to the 
lower IP on-off-pulse signal ratio. This can easily be seen in 
the off-pulse region of the average profile of the MP (Fig. 
lowest panel). 

The LRMI values increase by a factor of about four in 
the 1.4 GHz Lovell data sets and a factor of two in the 8.35 
GHz Effelsberg data sets for the two first components com- 
pared to those measured in session 1. Moreover, instead of 
the dip in the LRMI values shown in the previous session 
there is only a peak indicating high variability in the in- 
tensity of the observed pulses at the same pulse phase. The 
modulation indices for the Effelsberg IP tend to be lower to- 
wards the trailing edge of the IP profile. Due to the presence 
of the IP, there are two 2DFS plots, in Fig. H] The features 
visible in the vertically integrated 2DFS plots from session 
1 do not appear in session 2. 

In session 3 the 1.4 GHz Lovell data set (not shown) 
demonstrates a significantly large value of the LRMI at the 
trailing edge of the MP first component. This indicates in- 
frequent emission of very strong subpulses. The LRMI for 
the IP (not shown) was not calculated because of the large 
error in the measurement in modulation indices due to a 
small number of bright pulses. Because of the low S/N ratio 
in the 4.9 GHz WSRT observations, there are only a few val- 
ues of the LRMI calculated for the MP and none for the IP. 
The 8.35 GHz Effelsberg data set LRMI for the MP is two 
times lower than for the session 2 data set, hence is similar 
to the data set from session 1 (Fig. UJ right panel). The M2 
LRMI values become lower than those in Ml and M3. This 
indicates smaller intensity fluctuations of subpulses appear- 
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MP MJD 53886 MP MJD 53926 MP MJD 53933 




IP MJD 53926 IP MJD 53933 




Figure 4. Fluctuation analysis results shown from data obtained with the EfFelsberg radio telescope at frequency of 8.35 GHz during 
all 3 sessions. The results from different longitude ranges covering the main pulse and the interpulse are shown in the upper and lower 
rows respectively. For MP from session 2 (MJD 53926) only first two components are shown. The notation is the same as in Fig. [3] 



ing in the middle component and larger variations on the 
side components throughout the observation. The LRMI for 
the IP rises until it peaks in the middle of the average IP 
profile, which means that there are large intensity variations 
at this longitude. 



3.3 Pulse-energy distributions 

The results of the fluctuation analysis described in the pre- 
vious section clearly give evidence of intensity fluctuations 
over short periods of time, even within the same observing 
session. To study the pulse-amplitude characteristics of XTE 
J1810-197 we have made pulse-energy distributions for all 
frequencies and epochs for both the MP and the IP. Then, 
for each of the components of the MP and for the IP we fit 
their pulse-energy distributions with two well known models. 
Because the emission of XTE J1810-197 shows strong and 
spiky giant-pulse-like subpulses, we used p ower-law statis- 
tics to model this behaviour as shown by iLundgren et al.l 
l| 19951 ). while for the 'nor mal' pulses the logn ormal statistic 
was used as discussed bv lCairns et all l|200ll) : 



[E) oc , 
<E> 



lognormal 



2TTaE 



exp 



(1) 



(2) 



Since the power-law distribution integral is infinite, we 
introduced the parameter Emin, a minimum energy of the 
pulses, which together with p is used for fitting the power- 
law. The lognormal model is fit with /x and a. The average 
energy of the model distribution, < E > was set to match 
that of the observed distribution. In our analysis, we also 
took into account the effect of the noise. Since the noise in 
some of our observations was not pure Gaussian, due to the 
RFI present (e.g. 8.35 GHz Effelsberg data sets) we con- 
volved the observed noise with the model distributions. The 
derived parameters were optimised for the best fit by min- 
imising the using the Amoeba algorithm (|Press et al.l 
and are compared in Table (2] The detailed descrip- 
tion of the p roce dure of the fitting is w ell described by 
ICairnd (|2004l ) and lWeltevrede et al.l l|200^ . We discuss spe- 
cific cases below, but we note here that in general it was not 
possible to fit the distributions well with a single model. 
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Figure 5. Main pulse pulse-energy distributions from the 1.4 GHz Lovell (bottom row), 4.9 GHz WSRT {middle row) and 8.35 GHz 
Effelsberg {top row) observations from session 1. The components of the main pulse region are compared starting from first (Ml; left 
column) through second (M2; middle column) ending with third component (M3; right column). The horizontal axis denote energies 
which are normalised and scaled to the average pulse energy of the component. The off-pulse energy distributions of the components are 
shown with dashed lines. 



Fig. [S] presents the main pulse energy distributions for 
all the telescopes from session 1. The lowest row shows pulse- 
energy distributions from 1.4 GHz with ascending frequency 
towards the top row. The components of the main pulse re- 
gion are compared starting from Ml (left column) through 
M2 (middle column) ending with M3 (right column). For 
this session, only the Ml from the Lovell and WSRT data 
sets is best fit by lognormal pulse-energy distributions. The 
remaining components and the whole 8.35 GHz Effelsberg 
data set follows power-law-like statistics. The Effelsberg Ml 
is very strong and separates itself from the noise in the 
pulse-energy distribution (Fig. [S] top left-most plot). We 
emphasise that the Effelsberg Ml pulse-energy distribution 
is clearly a power-law without weaker underlying emission. 
This is different from the case of normal radio pulsar emis- 
sion, where the of giant pulses result in power-law tail ex- 
tending from the pulse-energy distribution. 

In session 2 the 1.4 GHz Lovell pulse-energy distribu- 
tion is best fit by lognormal behaviour. However, it also 
shows a long tail in the Ml pulse-energy distribution. Ml 
and M3 in the 8.35 GHz Effelsberg data set show similar 



behaviour (Fig. [6}. This introduces some complexity to the 
pulse-energy distribution fitting. Those pulse-energy compo- 
nents with long tails, cannot be fit with a single fit of either 
power-law or lognormal distribution. To try to get better 
fits, the distributions were split into two halves, the weak 
and strong pulses in terms of intensity, and then both halves 
of the distribution were fit with either lognormal or power- 
law distributions, like the giant pul se tails of normal radi o 
pulsar distributions are usually fit (|Lundgren et al.lll995l ). 
The Effelsberg IP (Fig. [T] left plot) shows a peak and a flat 
tail in its pulse-energy distribution, but it is possible to fit 
its pulse-energy distribution with a lognormal pulse-energy 
distribution. 

Session 3 is the only one, where the IP becomes present 
in all the data sets. Unfortunately, the IP subpulses from 
the 1.4 GHz Lovell data set are not seen above the noise 
level in the pulse-energy distribution, despite their appear- 
ance in the pulse stack. This is due to the weak S/N ratio 
and interference present in that data set. The Ml from the 
Lovell is best fit by a lognormal distribution, while M2 is best 
fit by a power-law distribution. The 4.9 GHz WSRT pulse- 
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Figure 6. Main pulse energy distributions from the 8.35 GHz Eflelsberg radio telescope observations from session 1 {bottom row), session 
2 (middle row) and session 3 {top row). The components of the main pulse region are compared starting from first (Ml; left column) 
through second (M2; middle column) ending with the third component (M3; right column). The horizontal axis denote energies which 
are normalised and scaled to the average pulse energy of the component. The off-pulse energy distributions of the components are shown 
with dashed lines. 




Figure 7. The interpulse energy distributions from the 8.35 GHz Effelsberg radio telescope observations from session 2 {left) and session 
3 {right). The off-pulse energy distributions are dashed lines. The horizontal axis denote energies which are normalised and scaled to the 
average pulse energy of the component. The off-pulse energy distributions of the components are dashed lines. 
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Table 2. The parameters of the best fits of the model distributions. Besides power-law, 
p and lognormal, /i, a distribution parameters, also the x'^i the number of degrees of 
freedom, ^.f a^nd significance probability, P{x^) are shown. 
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energy distributions show a power-law-like behaviour in all 
the components. A very strong IP is visible in the 8.35 GHz 
Efltelsberg data set, its pulse energy distribution is similar to 
the data set from session 2. The MP and the IP pulse-energy 
distributions show lognormal-like behaviour. Similarly to the 
previous session both the M2 and M3 pulses show a complex 
pulse energy distribution and the IP shows a peak, indicat- 
ing many pulses with flux close to the average energy of the 
IP, and a long tail in the pulse-energy distribution. Espe- 
cially the Effelsberg data set shows long tails. Considering 
only the Effelsberg strong pulses above 5 times the aver- 
age energy of the component, their pulse-energy distribu- 
tion is statistically insignificant because of the low number 
of pulses. The fit to the IP data from Effelsberg is also not 
significant. 

To better illustrate the strong intensity fiuctuations 
of XTE J1810-197 and to investigate any longitude- 
dependence of its strong subpulses, we have produced con- 
tour plots of the longitude-resolved cumulative pulse-energy 
distributions. Fig. |8] (left column) presents the distributions 
calculated for all the telescopes from session 1. The average 
profile is drawn with the thick line, while different linestyles 
denote the contours of the cumulative pulse-energy distribu- 



tion at the energy levels of 1%, 10% and 50%, normalised by 
the average peak energy < Ep > oi the average pulse profile. 
The 10% contour, for example, shows the energy level of the 
pulse-energy distribution above which 10% of pulses can be 
found. The second type of plot (right column; Fig. |8]) indi- 
cates the brightest time samples of each pulse-longitude bin 
(dashed line) scaled to the average intensity at that pulse- 
longitude bin, plotted on top of the average pulse profile 
(solid line). 

In the 1.4 GHz Lovell observation of session 1, sub- 
pulses appearing in the first component do not vary much, 
which indicates stable emission, this behaviour changes in 
M2 (Fig. [H] left column, top row) and this behaviour is re- 
peated in session 2 where it becomes very weak, although 
there are still some spiky subpulses present. The 4.9 GHz 
WSRT observation shows strong subpulses in Ml only, and 
the emission is stable throughout the rest of the profile dur- 
ing all sessions. 

To illustrate the session-to-session changes we used ob- 
servations from the Effelsberg telescope. The strong signal 
and the presence of the IP, as can be seen in Fig. (5] show 
how the emission from different components changes in time. 
Ml in the MP is always strong and the contour plots show a 
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Figure 8. The contour plot of the longitude-resolved cumulative pulse-energy distributions from the 1.4 GHz Lovell {top row), 4.9 GHz 
WSRT {middle row) and 8.35 GHz Effelsberg {bottom row) radio telescope observations from the session 1. Left: the thick solid line 
is the average pulse profile. The dashed line shows the brightest time sample for the each pulse-longitude bin. The other lines are the 

contours of the longitude-resolved cumulative pulse-cncrgy distributions at the energy level of the 1%, 10% and 50% compared to the 
average peak energy < Ep > of the average pulse profile. Right: the thick solid line is pulse profile. The dashed line shows the brightest 
time sample for each pulse-longitude bin compared with the average intensity < Ef,i„ > at that pulse longitude. 



strong stable emission (shown by 50% energy level contours) . 
It is also worth noticing that there are two peaks emerging 
in time. M2 grows stronger from session to session, and it 
always shows stable and strong emission. M3 is weak and is 
dominated by the presence of a few strong spiky subpulses. 



3.4 Single pulse correlation 

The last stage of our analysis was to perform a single pulse 
correlation. Before the analysis, we needed to align the data. 
The reason for this is the interstellar dispersion which causes 
the pulses received at higher frequencies to arrive earlier 
than the lower frequencies and the different path-lengths to 
the telescopes. Therefore all pulse arrival times have been 
corrected to a common reference frame. We transfered the 
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Figure 9. The contour plot of the longitude-resolved cumulative main pulse energy distributions from the 8.35 GHz Effelsberg radio 
telescope observations from session 1 {top row), session 2 {middle row) and session 3 {bottom row). The notation is the same as in Fig.|8] 



times of the pulses to the solar syste m barycentre, using 
the DE200 ephemeris (|Standish|[l98^ . Then, the overlap- 
ping (data sets at the different frequency pairs were taken 
for alignment, which was done by comparing the times of 
arrival of the pulses. Each data set was aligned in time with 
a one phase bin accuracy for each frequency pair. With the 
time resolution of 5.4 ms, the data sets were sufficiently 
aligned and all the telescope-specific delays were negligi- 
ble. We applied two techniques of single pulse correlation, 
the longitude-resolved linear correlation (LRLC) and the 
longitude-resolved cross-correlation (LRCC). Both methods 



are complementary and give interesting results as we will 
show below. 

The longitude-resolved linear c orrela tion is based on a 
method first introduced by IPopovI (Il986l ). In order to pro- 
duce the so-called 'linear correlation map' of the pulse-to- 
pulse variations, the l inear correlation ar r ay C j., as pre- 
sented in the work of iKarastergiou et al.l (|200lh needs to 
be calculated according to the following formula: 



n X af^i X ag 



.^[/.(fc)xg,(fc)-</. Xfl, >],(3) 



where fi(k) and gi{k) are time series of a bin i from two 
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distinctive frequencies with k being the pulse number, while 
af,i agj being their standard deviations respectively. 

The region of the highest linear correlation is expected 
to fall on the region around a diagonal which represents the 
time and spatial scale along which the subpulse intensities 
are linearly correlated. The LRLC method provides a good 
insight into linear dependent intensity variations of individ- 
ual pulses at different observed frequencies. However it does 
not take into account phase, shape information and non- 
linear dependencies between the components of the pulse 
from the different frequencies. The presence of low-intensity 
but persistent RFI in one of the data sets causes the appear- 
ance of a wave-like pattern in the linear correlation map as 
can be seen in the lower panel of Fig. IIUI 

An example of the total intensity linear correlation be- 
tween the 1.4 GHz LoveU and 4.9 GHz WSRT (upper plot) 
and the 1.4 GHz LoveU and 8.35 GHz Effelsberg (lower plot) 
data sets can be seen in Fig. (TO) The linear correlation of 
Ml is always strong between all frequency pairs. Most of the 
correlation regions in our results are very narrow. This con- 
firms the results from previous analysis steps, of strong and 
moderately stable emission from Ml at all frequencies and 
at all epochs. Surprisingly, the linear correlation of M2 is 
not as high despite the results showing it as the region with 
lowest intensity fluctuations which is characterised by low 
modulation indices in all our analysis. It is also interesting 
to note that the correlation seems to increase again towards 
the trailing edge of the pulse profile. The IP is visible only in 
the last two sessions, however in the case of session 2 the IP 
appears only in the 8.35 GHz Effelsberg data set. This con- 
strains the IP correlation analysis to the last session only. 
Unfortunately, the IP in the 1.4 GHz LoveU data set is very 
weak and the correlation results comprising this data set are 
unreliable. The linear correlation of the remaining data sets 
is equally strong across the whole IP range and higher than 
that of MP. 

Information on the shape and non-linear dependen- 
cies can be obtained by applying the LRCC method. This 
method of analysis is also more resistive to the presence of 
the periodic RFI in the 8.35 GHz Effelsberg data. The prin- 
ciple of producing the cross-correlation map is somewhat 
similar to the LRLC but instead of Eq. [3] we first apply the 
Fast Fourier Transform to both data sets. The next step is 
to multiply the first transformed data set by the complex 
conjugate of the second one. This produces a complex result 
of which the imaginary part is zero. Taking only the real 
part of the resulting cross-correlation coefficients matrix, 
the cross-correlation map is made as in the LRLC method. 
Examples of the results from the LRCC are presented in 
Fig. 111! The results of the cross-correlation analysis are sim- 
ilar to that of linear correlation. In session 3 we noticed the 
presence of a weak cross-correlation between the different 
components. Further analysis revealed that the occurrence 
of regions of higher correlation outside diagonal is caused 
by the RFI present in the 8.35 GHz Effelsberg data set. The 
region of correlation along the diagonal is very narrow in 
all of the correlated data sets. Similarly to the LRLC the 
cross-correlation is strong in Ml in all of the LRCC results, 
and can be explained by the similarities in shape and of the 
components at different frequencies. The next two compo- 
nents M2 and M3, are well correlated only between 4.9 and 
8.35 GHz data sets (Fig. 1111 right panel). The correlation 



results of 1.4 GHz with 4.9 or 8.35 GHz data sets do not 
show any correlation in M2 and M3 components. 



4 SUMMARY AND CONCLUSIONS 

We have analysed the data from the simultaneous multi- 
frequency single-pulse observations of XTE J1810-197 con- 
ducted at the frequencies of 1.4, 4.9 and 8.35 GHz, during 
observing sessions in May and July 2006. The phenomena 
revealed by our analysis indicate that the radio emission is 
clearly differe nt to the known radio pulsar propert i es. Pre - 
vious work bylCamilo et al.l ll2007bl'). ICamilo et al.l (|2007d '). 
iKramer et al.l (2007) and Lazaridis et al.l 1 20081 ) has shown 
an interesting overview of peculiarities of the pulsed ra- 
dio emission from the magnetar, e. g. flat spectral index, 
a — 0.0 ± 0.5, high degree of polarisation or long-term evo- 
lution of the polarisation angle swing with time. The results 
presented in our work confirm that the mechanism respon- 
sible for the magnetar radio emission appears to have a dif- 
ferent origin or perhaps even multiple origins, compared to 
the normal radio pulsars. 

XTE J1810-197 has a broad, multi-component profile 
with the IP only becomi ng present after session 1. Results 
from lKramer et al] (|2007l ) show that the separation between 
the MP and the IP is less than 180 deg which may be caused 
by an extremely wide MP beam (p ~ 44 deg) or by a non- 
dipolar magnetic field structure. The evolution of the aver- 
age pulse proflle takes place on day-to-day time-scales. This 
is clearly visible in the highest frequency data where the flrst 
component decreases in i ntensity, while the se cond and third 
become more prominent . iKramer et al.l l|2007l ) postulate, af- 
ter examining the average pulse profiles of XTE J1810-197 
from a greater number of observing sessions, that the mag- 
netar requires an unusually long observing time to obtain a 
stable profile or even lacks one, pe rhaps like PSR B0656-I- 14 
known for its bursting behaviour (jWeltevrede et al]|2006t ). 

The modulation index values change dramatically from 
one component to ano ther e ven w ithin a single observing 
frequency. iWeltevrede et al.l (|2007l ) define the modulation 
index to be the longitude-resolved modulation index rrii at 
the pulse longitude bin i where rrii has its minimum value. 
In most of the sources used in their analysis the LRMI shows 
a minimum in the middle of the pulse proflle with a typical 
value of m ~ 0.5, where the total intensity is relatively high. 
We also find that tendency in our results, but the magnitude 
of the values is somewhat closer to the values reported for 
the Crab Pulsar (m = 5). This agrees with the conclusion 
that due to the infrequent occurrence of the strong subpulses 
with narrow widths and broad distribution of the subpulse 
intensities, the modulation index is on average significantly 
larger in XTE J1810-197. We must note however, that in 
few cases, frequent occurrence of many strong and narrow 
subpulses with similar intensities at certain pulse phases of 
MP, results in lower LRMI values. The lowest values of the 
modulation indices occur during session 1, but even in this 
session they increase with increasing frequency as can be 
seen in Fig. |3] (upper panels). Such behaviour is in contrast 
to the normal radio pulsars where the LRMI values from 
lower freque ncies are on average hig her than that at higher 
frequencies (jWeltevrede et aLll2007l ). As we move in to ses- 
sion 2 the intensity fluctuations grow stronger with mini- 
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Figure 10. The plots of linear correlation between the 1.4 GHz Lovell and 4.9 GHz WSRT {left), and the 4.9 GHz WSRT and 8.35 
GHz Effelsberg data sets (right) from session 1. The contrast of the image was increased by showing only part of the cross-correlation 
coefRcient range. The bottom and left panel hold the average profiles of the correlated data sets. The bar on the right side of the plot 
denotes the different levels of the linear correlation coefficient. The value of -0.2 means 20% negative linear correlation while 0.85% 
positive linear correlation is denoted by 0.85 value. The maximal and minimal values of the linear correlation coefficients are shown in 
the left lower part of plots. Characteristic wave-like pattern visible in the lower plot is a result of the correlation of the signal from the 
Lovell with the RFI present in the Effelsberg data set. 
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Figure 11. The plots of cross-correlation between the 1.4 GHz Lovell and 8.35 GHz Effelsberg (left), and the 4.9 GHz WSRT and 8.35 
GHz Effelsberg data sets (right) from session 3. The contrast of the image was increased by showing only part of the cross-correlation 
coefficient range. The profiles in the bottom and left panel represent the average profiles. The bar on the right side of the plot denotes the 
different levels of the cross-correlation coefficient normalised for clarity. The value of 0.0 means no cross-correlation while 50% positive 
cross-correlation is denoted by 0.5 value. The maximal and minimal values of the cross-correlation coefficients are shown in the left lower 
part of the plots. 



mum modulation indices at values of around 4 for the 1.4 
GHz Lovell data sets. For the 8.35 GHz Effelsberg data, we 
calculate the minimum modulation indices to be 1 in the MP 
and 2 in the IP region. It is remarkable, that for both data 
sets in this session the modulation indices in the MP peak 
at values of 7.5 and 10 for the Lovell and Effelsberg data 



respectively. These values are extremely high and except for 
the Crab Pulsar, unprecedented in the re sults from modu- 
lation analysis from normal radio pulsars (jWeltevrede et al] 
1200 it ). In session 3, for the 4.9 GHz WSRT data sets due to 
the low signal-to-noise ratio, the LRMI values are not sam- 
pled densly throughout the whole pulse profile range, but 
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available values are comparable with session 1. The Lovell 
and EflFelsberg LRMI values are similar to that of session 
1. The longitude-resolved modulation analysis results pre- 
sented above show the variation of the LRMI values on day- 
to-day time-scales and dramatic changes with pulse phase. 
In all session we find the frequency evolution of the LRMI 
values in contradiction with proprieties of normal radio pul- 
sars. 

The following step of our analysis, the 2DFS, do not 
show any regular drifting behaviour in any sessions or any 
frequencies. However, we find phenomenon manifesting itself 
as characteristic bumps in the XTE J1810-197 vertically col- 
lapsed 2DFS. We interpret this phenomenon as the tendency 
of the subpulses to be equally separated in the consecutive 
pulses' profiles throughout observation. We also notice the 
presence of a peak visible in the vertically collapsed 2DFS at 
1.4 GHz LoveU, 4.9 GHz WSRT from session 1 (Fig. H left 
and middle panels) and 8.35 GHz Effelsberg from session 3 
(Fig. |4j right panel in the upper row) data sets. We interpret 
the peak as a signature of the baseline variations in those 
data sets. The lack of regular drift from the magnetar might 
be associated with its rapidly changing e mission properties 
and yo ung age (r < lOkyr). In their work. IWeltevrede et al.l 
(12003) have shown that the fraction of young pulsars show- 
ing regular drifting is very low. Although, one could also ar- 
gue that the strong radio variability might mask any regular 
structures or the physics of the magnetar's radio emission is 
different from radio pulsars. 

The high variability of the magnetar emission is also 
reflected in its pulse-energy distributions. For all sessions at 
all observed frequencies we have made pulse-energy distri- 
butions for each of the main-pulse components as well as for 
the IP (whenever present). We fit each of the pulse energy 
distributions with models based on a power-law or lognor- 
mal statistics for comparison between our observations and 
existing pulsar emission models. As justified later in this 
section, propagation effects in the interstellar medium are 
negligible for our data analysis and are therefore ignored. 
Table [2] presents the results of the best fits. The significance 
of the best fits are low due to the oversimplified models, but 
changes in the best-fit models of the components in different 
sessions are nevertheless, very peculiar. We interpret that as 
indicating the possible presence of multiple emission mech- 
anisms with different statistical behaviour embedded in the 
same pulse phases. 

In a series of papers. Cairns and collaborators 



JCairns et al.[l200ll . ICairns et al1l2003al . ICairns et alll2003bl . 
ICairns et all 120041 1 investigate different possible models of 



the emission physics using observations of PSRs B0833-45 
(Vela Pulsar), B1641-45 and B0950-I-08. They show that 
the models of emission vary with changing pulse phase 
in the analysed sources. They find in these pulsars, that 
longitude-resolved pulse-energy distributions which appear 
to be lognormal and are representative of the normal pul- 
sar radio emission, can be fit with a single emission model, 
or convolution of Gaussian-lognormal or double lognormal 
models. This multi-model manifestation is explained as two 
non-related waves coupling together in the inhomogeneous 
plasma in the pulsar magnetosphere. While Carins et al. 
successfully fit the longitude-resolved pulse-energy distri- 
butions with lognormal statistics, in many cases it is also 
valid to use it to fit the integrated pulse-energy distributions 



(iJohnston fc Romani|[2003 ). The presence of approximately 
power-law distributions is caused by these different phenom- 
ena. The origin of these phenomena can be identified based 
on the values of power-law indices and can be associated 
with (i) giant pulses i. e. pulses which have integrated fiux 
density greater than 10 times the integrated fiux density 
of the average profile, (ii) giant micropulses, which have a 
peak flux density that is 10 times the peak flux density of 
the average profile at the same phase, but less than 10 times 
integrated flux density of the average profile, (iii) so-called 
precurs or emissions occurring just before the MP as for Crab 
Pulsar (jCordes et al . 2004, M offett. D. A. fc Hankins. T. H.I 
Il996l ). which have power-law distributions that are non in- 
trinsic and caused by Gaussian distributions from the back- 
ground normal pulsar emission convolved with higher en- 
ergy lognormal components. We argue that in the case of 
the magnetar. Ml consists of strong and narrow giant-like 
subpulses followed by a components of weak precursor-like 
subpulses. 

The giant pulse and giant micropulse phenomena are 
known in a few rad io pulsars like PSRs B053H-21 (Crab 
Pulsar: IStae lin fc Reifcnstein 1968 1. B0833-45 (Vela Pulsar ; 
I Johnston et~ l200ll). B1937-H21 (|Soglasnov et all 120041 ) 
or B1133-I-16 ^Kramer et al.ll2003l ). The emission of these 
pulses can be characterised by the power-law energy distri- 
butions, broadband emission or occurrence within narrow 
pulse phases. These phenomena are believed to be associ- 
ated wi th the high-energy emission in the out er magneto- 
sphere (|johnston fc Romani|[20o3 . ICairnsll2004h . In the case 
of XTE J1810-197 we also see strong spiky subpulses which 
could be associated with the giant pulse phenomenon, but 
their widths are larger than that of giant pulses of normal 
radio pulsars. Also, their occurrence, which covers the whole 
longitude range of that component, stand in contradiction 
to this deflnition. The most prominent example illustrat- 
ing the above phenomenon in our observations is the first 
component, especially in the Effelsberg data sets from ses- 
sions 2 and 3 as can be seen in Fig. [6] This component has 
many strong and spiky subpulses appearing within its whole 
longitude range, which dominates the high-energy tail in 
its pulse-energy distributions. This makes fitting the pulse- 
energy distributions with only one law impossible. However, 
the attempt to decouple the component distribution into 
low and high energy parts also did not result in good fits. 
A sim ilar case occurs for PSR B0656-I-14 (Weltevrede et akl 
[2003), where weak emission is coupled with a component re- 
sponsible for high energy bursts. In the case of XTE J1810- 
197 there may be more than two models contributing at 
one phase in the observed distributions. As shown later the 
magnetars' very dynamic magnetosphere may be an expla- 
nation for such multi-component pulse-energy distributions. 
This argues that the emission is in general broadband, but 
the degree of variability is very different. Those components 
average together in pulse-energy distributions, which makes 
them difficult to fit properly using known statistical models. 

Despite the changes in the pulse profile and pulse- 
amplitude characteristics on short time-scales, the correla- 
tion analysis gave results which contradict the overall pic- 
ture of unstable emission from XTE J1810-197. The LRLC 
analysis shows significant correlation results in the majority 
of the frequency pairs used. The narrow and high correlation 
regions denote significant dependence between the intensi- 
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ties of the subpulses on small time- and spatial scales. The 
correlation always occurs between the first components in 
all of the analysed frequency pairs, with sporadic correlation 
between the third components. In contrast, the second com- 
ponent, is found to be a stable emission region with lower 
modulation indices, was very weakly correlated. To examine 
the non-linear dependence between frequencies we used the 
LRCC method. The correlation is weaker when compared 
to the LRLC method, but correlated regions are also very 
narrow, showing that there are similarities in the phase and 
shape of the subpulses at different observed frequencies. 

While the larger time-scale flux variations might be 
normally attributed to interste llar scintillation, they have 
been rejected i n the work of ICamilo et al.l (|2007d ) and 
iLazaridis et al.l l|2008t ) as being responsible for the variability 
over short time scales. This points to the intrinsic behaviour 
of the magnetar as a cause. The lack of a regular drift, broad 
pulses, the presence of subpulses with quasi-periodic mod- 
ulation, difficulties with fitting the data with single lognor- 
mal or power-law models allow us to draw a conclusion of 
non-stable emission due to the possible turbulent magnetar 
magnetosphere. A model explaining that emissi on has been 
propo sed very recently. In his work, Thompson l|Thompsonl 
l2008al lb[) gives an extensive explanation of the pair creation 
processes in ultra-strong magnetic field and particle heating 
in a dynamic magnetosphere. He considers the details of the 
QED processes that create electron-positron pairs in high 
magnetic fields of the order of 10^'' G. He discusses the pos- 
sibility of a strong enhancement of the pair creation rate in 
the open-field circuit and outer magnetosphere by instabili- 
ties near the light cylinder. Thompson also refers to the flat 
radio spectra as a possible result of the high plasma density 
in the open magnetic held lines. One of the model expla- 
nations of the magnetar's broad pulse profile, is its beam 
geometry. In normal radio pulsars, wide pulse profiles are 
usually caused by the alignment between rotation and mag- 
netic pole axis. The line of sight of the observer stays within 
the emission beam for a large fraction of the pulse period 
resulting in the long duty cycle. In the case of XTE J1810- 
197 the solution of fittin g the position angle swing with the 
Rotating Vector Model jRadhakrishnan fc Cookg|l969l ) re- 
sults in non-aligned geometry (a = 44 deg, /3 = 39 deg), but 
the beam radius inferred fro m the MP pulse has a width of 
about p ~ 44 deg as shown bv lKramer et al.ll|2007l ). This re- 
sult excludes viewing geometry as a reason for a wide pulse 
profile in XTE J1810-197. The model of the dynamic outer 
magnetosphere has a promising application in explaining the 
radio emission from the magnetars and is consistent with 
the magnetars' emission features such as flat radio spectra, 
broad pulses and rapid variability. 

Since the detection of radio emission from XTE J1810- 
197, its relation to the new clas s of objects called Ro - 
tating RAdio Transients (RRATs. [McLaughlin et al.l 120061 ) 
was suggested. The lack of X-ray counterpart in any of 
kn own RRAT sources a rgued against this hypothesis un- 
til [Reynolds et all l|2006l ) reported the flrst X-ray detection 
from RRAT J1819-1458. This RRAT was known to emit ra- 
dio bursts with 4.26 s spin period based on its timing anal- 
ysis and the inferred dipole surface magnetic field of about 
5 X 10^^ G. The detection of t he periodic X-ray pulsati ons 
aligned with the radio bursts bv lMcLaughlin et al. I i2003) al- 
lowed the comparison of its X-ray emission proprieties with 



XTE J1810-197 and excluded a close relationship with this 
RRAT. However , it would be very interesting as proposed by 
iRea et al.1 (|2008l ) to search for the RRAT-like radio bursts 
from XTE J1810-197 at its quiescence level. This could in- 
vestigate the hypothesis of a link between the magnetars, 
RRATs and young radio pulsars. 
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